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INTRODUCTION 

This  final  report  on  Contract  N0014-74-C-0273  summarizes  work  on 
the  bottom  features  task  during  the  four  year  period  of  this  contract, 
concentrating  especially  on  the  fourth  year  results.  A report  on  the 
beach  environment  task  of  this  contract  has  been  prepared  under  separate 
cover.  The  work  was  accomplished  at  the  Environmental  Research  Institute 
of  Michigan  under  the  guidance  of  Principal  Investigators  Robert  Vincent 
and  later  Fred  Thomson.  The  performance  period  for  this  contract  was 
January  1974  - February  1978. 

This  introduction  contains  a description  of  the  bottom  features  task 
and  a discussion  of  its  relevance  to  the  coastal  reconnaissance  problem. 

A summary  of  the  accomplishments  made  during  this  contract  is  contained 
in  Section  2,  and  a detailed  technical  discussion  follows  in  Section  3. 

The  bottom  features  task  arose  from  earlier  research  on  water  depth 
mapping  using  the  selective  penetration  of  light  in  water  in  various 
spectral  bands  as  sensed  by  a passive  multispectral  scanner.  Since  the 
reflected  light  is  influenced  by  both  the  water  depth  and  the  bottom 
type,  these  two  problems  are  closely  related.  In  the  case  of  depth 
mapping,  the  object  is  to  combine  the  signals  in  various  wavelength 
bands  in  order  to  produce  an  output  signal  which  varies  only  with  depth, 
independently  of  the  bottom  type.  Conversely,  in  the  case  of  bottom 
features  mapping  the  problem  is  to  remove  the  effect  of  depth  variations 
so  as  to  produce  a depth-invariant  indicator  of  bottom  type.  It  now 
appears  that  the  bottom  features  problem  is  the  prior  one  in  the  sense 
that  bottom  type  information  can  be  extracted  without  knowledge  of  depth, 
but  depth  information  cannot  be  reliably  extracted  in  the  general  case 
without  knowledge  of  the  bottom  type. 

Research  on  the  extraction  of  water  depth  information  from  multi- 
spectral scanner  data  has  been  carried  out  at  ERIM  since  1967  under  the 
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sponsorship  of  of  various  federal  agencies  including  NASA,  NOAA/NESS, 
and  the  Defense  Mapping  Agency.  Techniques  for  extracting  bottom-type 
information  have  also  been  studied  at  ERIM  since  1972,  first  under  the 
sponsorship  of  the  Environmental  Protection  Agency  and  since  1974  for 
the  Office  of  Naval  Research.  The  specific  purpose  of  the  research 
sponsored  by  EPA  was  to  map  the  distribution  of  Cladophora  in  Lake 
Ontario,  while  the  efforts  supported  by  ONR  have  beey  directed  toward 
the  general  problem  of  mapping  bottom  features. 

The  relevance  of  the  bottom  features  mapping  task  to  the  general 
coastal  reconnaissance  problem  lies  not  only  in  its  connection  with 
depth  mapping,  but  also  in  the  information  it  can  give  about  coastal 
geology  and  ecology.  Information  about  the  subsurface  geology  can  some- 
times be  given  directly  (e.g.,  through  differentiation  of  sand,  rock, 
mud,  coral,  etc.)  and  sometimes  indirectly  through  the  association  of 
benthic  algal  communities  with  different  substrate  types.  Information 
about  coral  or  algal  communities  is  useful  in  ecological  studies, 
relating  to  natural  marine  ecosystems  and  the  impact  of  waste  discharge, 
runoff,  and  resource  extraction  activities. 


2 


2 


SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  first  work  done  at  ERIM  on  the  problem  of  bottom-features 
recognition  was  based  on  an  analysis  of  imagery  generated  by  an  analog 
ratio  processing  technique,  using  multispect ral  data  collected  by  ERIM 
aircraft  over  Lake  Ontario  and  the  east  coast  of  Florida.  The  disadvantage 
of  this  technique  was  that  the  processing  was  not  able  to  be  precisely 
controlled  and,  therefore,  was  not  repeatable.  During  the  first  year 
of  the  ONR  program  (1974)  an  analysis  was  carried  out,  using  a simple 
water  reflectance  model,  which  resulted  in  the  definition  of  a modified 
ratio  algorithm  (MRA)  and  the  implementation  of  this  algorithm  on  a 
digital  computer.  During  tlve  second  year  (1975)  an  investigation  of 
the  limitations  of  the  MRA  was  begun  using  the  simple  radiance  model, 
and  a more  accurate  model  incorporating  the  effects  of  scattering  in 
the  water  was  developed.  In  the  third  year  (1976)  the  model  was  extended 
to  include  atmospheric  effects,  and  a more  thorough  evaluation  of  the 
MRA  was  accomplished  using  this  model.  Tills  analysis  indicated  certain 
basic  limitations  of  the  MRA  and  suggested  a more  general  algorithm 
based  on  a logarithmic  transformation  of  the  radiance  data  collected 
by  the  scanner. 

The  limitations  of  the  modified  ratio  algorithm  are  that  (1)  it 
is  inherently  a two-band  technique,  (2)  the  wavelength  bands  must  be 
chosen  carefully  for  each  water  type,  such  that  the  water  attenuation 
coefficients  are  the  same  in  both  bands,  and  (3)  the  technique  cannot 
discriminate  between  hot  tom  types  with  equal  reflectance  ratios  in  these 
bands.  As  a result  of  these  limitations,  the  average  probability  of 
correct  classification  for  three  typical  bottom  types  (sand,  mud,  and 
vegetation)  in  clear  coastal  water  (Jerlov  Type  3)  is  only  about  83% 
at  zero  water  depth  and  falls  to  50%  at  a depth  of  6 meters  for  a typical 
aircraft  scanner. 
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During  the  last  year  of  this  contract  an  Improved  bottom  recognition 
algorithm  was  developed  which  is  not  subject  to  those  limitations.  This 
algorithm  can  utilize  any  number  of  bands,  and  the  choice  of  these  bands 
is  not  limited  by  the  constraints  imposed  in  the  MRA.  In  fact,  by 
choosing  bands  with  different  water  attenuation  coefficients,  materials 
with  similarly  shaped  reflectance  spectra  which  could  not  be  separated 
by  the  MRA  can  be  discriminated  using  the  new  technique.  The  two-band 
version  of  this  algorithm  lias  undergone  limited  theoretical  and  empirical 
evaluations , which  are  the  subject  of  this  report.  For  the  situation 
described  above,  the  improvement  in  classification  accuracy  is  shown  in 
Figure  1. 

A theoret leal  water  reflectance  model  developed  for  the  purpose  of 
evaluating  the  bottom  recognition  algorithms  lias  proved  to  be  a valuable 
product  of  the  ONR  program  in  its  own  right.  This  model  has  been 
described  in  previous  reports  [1-3]  and  will  also  be  reviewed  in  Section 
3 of  this  report.  The  model  has  been  used  in  various  programs  for  the 
Navy  [4)  and  the  National  Oceanic  and  Atmospheric  Administration  [5-7], 
and  is  planned  for  use  in  several  new  programs  for  NASA.  The  results 
accomplished  during  this  contract  may  be  summarized  as  follows: 

1)  Through  a sequence  of  modeling  and  algorithm  development  steps, 
a multichannel  algorithm  has  been  developed  for  bottom  type  recognition, 
under  a variable  depth  of  water.  The  optimum  channels  for  this  algorithm 
are  selected  primarily  on  the  basis  of  spectral  differences  in  bottom 
type.  The  probability  of  correct  classification  for  a typical  situ- 
ation using  the  latest  algorithm  and  the  previous  algorithm  (MRA)  are 
shown  in  figure  1. 

2)  As  an  outgrowth  of  algorithm  development,  a radiative  transfer 
model  has  been  developed  which  is  able  to  calculate  the  spectral  radiance 
received  by  a sensor  at  any  altitude,  under  arbitrary  solar  illumination, 
view  angles,  and  atmospheric  conditions,  and  for  arbitrary  water  properties 
and  bottom  reflectances. 
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FIGURE  1.  AVERAGE  PROBABILITY  OF  CORRECT  CLASSIFICATION  OF  THREE  BOTTOM 
TYPES  (SAND,  MUD,  AND  VEGETATION)  IN  JERLOV  WATER  TYPE  3 USING 
NEW  BOTTOM  RECOGNITION  ALGORITHM  (UPPER  CURVE)  AND  THE  MOD- 
IFIED RATIO  ALGORITHM  (LOWER  CURVE). 
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3)  A comprehensive  evaluation  of  the  accuracy  of  the  new  algorithm 
requires  knowledge  of  coastal  water  absorption,  scattering,  and  scattering 
phase  functions  which  is  not  presently  available.  Therefore,  in  evaluating 
the  algorithms,  several  assumptions  about  these  properties  were  made. 

For  a more  comprehensive  evaluation,  in  particular,  for  evaluation  of 
candidate  coastal  sensor  designs,  a program  to  obtain  the  required 
information  should  be  initiated. 

These  conclusions  lead  to  the  following  recommendations  for  further 
study: 

1)  A quantitative  empirical  evaluation  of  the  bottom  recognition 
algorithm  should  be  made,  by  comparing  results  of  processing  multispectral 
scanner  data  over  Panama  City  and  the  Bahamas  with  actual  field  conditions. 

2)  A broader  base  of  information  on  important  water  and  bottom  parameters 
(i.e. , water  absorption  and  scattering  coefficients,  scattering  functions, 
and  bottom  reflectances)  should  be  acquired,  to  permit  more  realistic 
simulations  and  a more  complete  theoretical  evaluation  of  the  algorithm. 

3)  Theoretical  and  empirical  evaluations  of  the  multichannel  algorithm 
with  more  than  two  channels  should  be  carried  out.  For  a closer  explanation 
of  this,  refer  to  page  13,  paragraph  2,  of  this  report. 

4)  An  evaluation  should  be  made  of  the  applicability  of  this  algorithm 
to  satellite  data,  and  the  utility  of  this  application  should  be  examined 

in  view  of  the  spectral  channels  and  spatial  resolution  available  with 
current  satellites. 
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ALGORITHM  DEVELOPMENT  AND  EVALUATION 

This  section  contains  a discussion  of  the  bottom  recognition 
algorithm  developed  during  the  fourth  year  of  this  contract,  including 
the  methods  used  for  developing  and  evaluating  the  algorithm  and  the 
results  of  this  evaluation.  The  theoretical  water  radiance  model  used 
in  the  development  and  evaluation  of  the  algorithm  is  described  in 
Section  3.1.  The  algorithm  itself  Is  described  in  Section  3.2,  and 
the  evaluation  of  the  algorithm  theoretically  and  empirically  is  discussed 
in  Sections  3.3  and  3.4,  respectively. 

3.1  WATER  RADIANCE  MODEL 

The  water  radiance  model  used  for  developing  and  evaluating  the 
bottom  recognition  algorithm  is  a modification  of  the  quasi-single- 
scattering approximation  [8]  which  includes  the  effects  of  reflection 
from  the  bottom  and  from  the  air-water  interface,  as  well  as  scattering 
and  absorption  in  the  water.  This  approximation  has  the  advantage  of 
computational  speed  over  exact  numerical  methods,  and  gives  results 
which  are  probably  within  the  range  of  measurement  accuracy.  Comparisons 
with  exact  calculations  were  presented  in  an  earlier  report  [11,  and  the 
results  of  an  empirical  test  follow  in  this  section. 

An  atmospheric  model  [91  was  also  combined  with  the  water  radiance 
model  in  order  to  fully  simulate  the  radiance  received  at  an  airborne 
sensor.  This  atmospheric  model  is  used  to  calculate  the  solar  irradiance 
at  the  surface,  the  sky  radiance,  the  path  radiance,  and  the  atmosphor 1- 
transmittance.  The  total  radiance  at  the  sensor  aperture  Is  given  bv 

L - jn2T  L(li’  ,<j>)  + r lJ  T + L (3.1) 

tot  f s ssjap 
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where  n - index  of  refraction  of  water 

T = transmittance  of  air-water  interface 
s 

L(u',<t>)  » upwelling  underwater  radiance  (c.f.  Appendix,  p.  33) 

r = reflectance  of  air-water  interface 
s 

= sky  radiance  (from  atmospheric  model) 

T = atmospheiic  transmittance  (from  atmospheric  model) 
a 

Lp  = path  radiance  (from  atmospheric  model) 

An  experiment  conducted  jointly  by  ERIM  and  the  Naval  Coastal  Systems 

Center  (NCSC)*  in  May  1977  [10]  presented  an  opportunity  for  empirically 
evaluating  the  water /atmospheric  model.  In  s itu  measurements  were  made 
of  the  reflectance  of  three  bottom  types  in  St.  Andrew  Bay,  Florida,  using  an 
ISCO  spectroradiometer  and  three  known  reflectance  panels.  Samples 
of  bright  sand  from  the  Gulf  of  Mexico  beach  near  Panama  City,  Florida 
were  also  taken  and  laboratory  measurements  were  made  of  the  reflectance 
of  those  samples  at  ERIM  using  a Cary-14  spectrometer.  The  measured 
reflectances  for  these  four  bottom  types  are  shown  in  Figure  2.  Water 
optical  parameters  were  measured  ijn  situ  by  NCSC  personnel,  and  deter- 
mined independently  by  an  analysis  of  the  radiances  measured  by  the 
ERIM  M-8  airborne  multispectral  scanner.  Diffuse  attenuation  coefficients 
for  the  Gulf  and  Bay  waters  determined  by  the  latter  procedure  are  shown 
in  Figure  3.  These  values  agree  quite  closely  with  the  photometer 
measurements  (K  * 0.15-0.25  m ^ in  Gulf,  and  K = 0.35-0.45  m ^ in  Bay) 
made  by  NCSC  [11].  Additional  measurements  of  the  beam  attenuation 
coefficient  (a)  allowed  an  estimate  of  the  total  scattering  coefficient 
(s)  using  the  approximate  relationship 

s = a - K (3.2) 

This  estimate  of  s was  0.5  m * for  the  Gulf  and  1.0  m 1 for  the  Bay. 

Since  the  Bay  presents  the  more  interesting  situation  in  terms  of 
scattering  and  diverse  bottom  conditions,  a set  of  radiances  was 

* 
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calculated  for  the  Bay  using  the  parameters  discussed  above,  and 
assuming  an  average  particle  scattering  function  from  the  measurements 
of  Petzold  [12].  The  radiances  were  calculated  for  a water  depth  of 
1 meter,  and  compared  with  radiances  measured  by  the  M-8  scanner  at 
approximately  the  same  depth.  The  comparison  is  shown  in  Figure  4. 

The  only  systematic  difference  between  the  calculated  and  measured 
values  is  in  the  blue  region  ( .475- . 50pm) . It  is  not  known  whether 
this  difference  is  due  to  errors  in  measurement  of  parameters,  deficiencies 
in  the  model,  or  miscalibration  of  the  scanner.  However,  the  overall 
error  is  quite  small  and  the  comparison  was  considered  to  be  a confirma- 
tion of  the  model. 

3.2  BOTTOM  RECOGNITION  ALGORITHM 

The  bottom  recognition  algorithm  is  a technique  for  transforming  a 
set  of  N variables  (essentially  the  radiances  in  N wavelength  bands) 
which  contain  both  water  depth  and  bottom  type  information  Into  a set 
of  N-l  depth-invariant  variables  containing  only  bottom  type  information, 
and  one  variable  containing  depth  information.  This  process  is  actually 
a composite  of  two  separate  operations.  The  first  operation  is  performed 
by  subtracting  the  deep-water  signal  from  each  data  point  and  taking 
the  natural  logarithm  of  the  difference.  The  purpose  of  this  step  is 
to  linearize  the  depth-dependence  of  the  bottom-reflected  signals.  The 
water  radiance  model  described  in  Section  3.1  shows  that  this  lineariza- 
tion is  only  approximate,  but  that  significant  nonlinearities  only 
occur  for  very  bright  bottoms  and  shallow  depths.  The  second  operation 
is  a linear  transformation  which  has  the  properties  of  a coordinate 
system  rotation.  The  purpose  of  this  transformation  is  to  remove  the 
depth  dependence  from  all  but  one  of  the  signal  channels. 

The  mathematical  formulation  of  this  algorithm  is  detailed  in  the 
Appendix  of  this  report  for  the  general  case  of  N channels.  A method 
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Grass 


FIGURE  4.  CALCULATED  (CURVES)  AND  MEASURED  (POINTS)  RADIANCES  FOR 
WATER  DEPTH  ' 1 METER,  ST.  ANDREW  BAY,  FLORIDA. 


of  evaluating  the  performance  of  the  algorithm  for  the  two-channel  case 
is  also  described  in  the  Appendix,  and  results  are  presented  for  one 
specific  case.  Further  results  are  presented  in  Section  3.3.,  including 
a determination  of  optimum  bands  for  three  different  water  and  bottom 
types. 

For  the  two-band  case  the  algorithm  results  in  a single  depth-invariant 
variable  which  may  be  used  as  an  index  of  the  bottom  type.  For  the  case 
of  three  or  more  input  bands,  the  algorithm  results  in  two  or  more  depth- 
invariant  channels.  Thus,  a further  step  is  required  to  combine  this 
information  into  a single  index  of  bottom  type.  Such  multivariate  pat- 
tern recognition  techniques  exist,  but  they  have  not  as  yet  been  combined 
with  the  algorithms  described  in  this  report.  Consequently  only  the  two- 
band  case  has  been  fully  implemented  and  evaluated.  It  is  expected  that  the 
addition  of  a third  band  will  yield  a small  but  possibly  significant 
improvement  in  classification  accuracy,  and  that  the  incremental  improvements 
for  more  than  three  bands  will  decrease  with  the  number  of  bands  employed. 

3.3  THEORETICAL  EVALUATION 

A computer  program  was  written  to  perform  an  error  analysis  for  the  two- 
channel  bottom  recognition  algorithm,  using  the  radiance  model  described  above 
to  simulate  the  daca  collected  by  a multispectral  scanner  system.  The  program 
considers  all  possible  pairs  of  the  wavelength  bands  entered  and  computes  the 
probability  of  misclassi f icat ion  for  each  bottom  type  at  each  depth  entered. 

For  each  wavelength  pair,  the  program  calculates  the  coordinate  system 
rotation  parameters  as  described  in  the  Appendix  using  the  known  water  attenu- 
ation coefficients.  The  values  of  Y^  (the  depth-invariant  variable)  are  cal- 
culated at  the  "training  depth"  for  each  bottom  type,  and  the  decision  boun- 
daries are  chosen  midway  between  the  neighboring  values  of  Yj. 

After  the  "training"  process  is  completed,  the  radiance  data  is  read  again 
and  the  probability  of  classification  in  each  category  is  calculated  for  the 
value  of  NEAL  (the  noise-equivalent  radiance)  entered.  The  results  for  sand,  mud. 
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and  vegetation  in  Jerlov's  [13]  water  type  3 are  reported  in  the  Appendix 

-2  -1  -1 

for  a noise  figure  NEAL  - .05  mW  cm  sr  um  . This  value  is  typical  for 

an  aircraft  scanner  such  as  the  ERIM  M-7  and  M-8  scanners,  but  may  be 

improved  by  techniques  such  as  spatial  filtering  or  by  changes  in 

system  design.  In  the  1977  NCSC  experiment  [14]  the  noise  was  reduced 

-2  -1  -1 

to  approximately  .02  mW  cm  sr  wn  by  applying  a 3 x 3 pixel 
smoothing  function  to  the  data  collected  by  the  M-8  scanner. 


The  error  analysis  was  repeated  for  each  possible  wavelength  pair 
in  the  range  0.4-0.6*m  with  0.025um  intervals,  using  the  lower  noise 
figure.  The  results  for  water  type  3 are  shown  in  tables  1 and  2 for 
depths  of  2 meters  and  4 meters,  respectively.  These  results  were  obtained 
using  a training  depth  of  2 meters. 

The  results  shown  in  tables  1 and  2 reveal  two  distinct  regions 
where  the  performance  is  optimized.  The  band  pair  .475-525vim  is  repre- 
sentative of  the  first  region  and  the  band  pair  .525-.60um  is  representa- 
tive of  the  second  region.  Although  the  second  band  pair  gives  slightly 
better  results  at  2 meters,  the  first  band  pair  has  the  advantage  of 
being  less  sensitive  to  changes  in  water  quality.  Figure  5 shows  the 
average  percent  correct  classification  for  the  band  pair  .475-525um  as 
a function  of  depth  in  water  type  3,  assuming  the  algorithm  has  been 
trained  in  this  water  type  at  2 meters  depth,  and  the  performance  in 
Jerlov's  water  type  5 with  the  same  operating  parameters.  Thus,  the 
algorithm  correctly  classifies  the  bottom  about  80%  of  the  time  in  water 
type  3 at  4 meters  depth,  but  if  a portion  of  the  scene  contains  water 
type  5,  the  classification  accuracy  for  this  portion  drops  to  about 
50%  at  a depth  of  4 meters. 
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AVERAGE  CLASS  IE' l GAT  ION  ACCURACY  FOR  SAN1),  MUD,  AND 
VEGETATION  IN  WATER  TYPE  3 Al  2 METERS  DEPTH 
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TABLE  2 

AVERAGE  CLASSIFICATION  ACCURACY  FOR  SAND,  MUD,  AND 


VEGETATION 

IN  WATER  TYPE 

3 AT  4 

METERS  DEPTH 
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38 

33 
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.525 
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Percent  Correct  Classification 


Similar  calculations  were  carried  out  for  water  type  5,  training  the 
algorithm  in  this  water  type  at  2 meters  depth.  The  results  of  these 
calculations  are  shown  in  tables  3 and  4.  The  optimum  wavelength  pairs 
seem  to  be  very  nearly  the  same  for  this  case,  suggesting  that  the 
optimum  wavelengths  are  determined  primarily  by  the  bottom  types  rather 
than  by  the  water  type.  The  decision  boundaries  selected  for  water 
type  5 are  similar  to  those  for  water  type  3 in  the  first  group  (around 
.475-.  525vjm)  but  are  quite  dissimilar  in  the  second  group  (around  .525-. 60pm). 
This  accounts  for  the  relative  insensitivity  of  the  algorithm  to  changes 
in  water  quality  using  the  first  wavelength  pair,  and  the  large  sensiti- 
vity using  the  second  wavelength  pair.  Figure  6 shows  the  performance 
of  the  algorithm  trained  in  water  type  5,  using  the  band  pair  .475-. 525pm, 
and  evaluated  for  water  type  5 and  water  type  3. 

A simulation  was  also  done  for  the  St.  Andrew  Bay  test  site,  with 
sand,  shoal  grass,  and  turtle  grass  as  the  bottom  types.  The  water 
attenuation  coefficients  at  this  site  (c.F.  Figure  3)  were  intermediate 
between  Jerlov  types  5 and  7.  The  algorithm  was  trained  at  1 meter  and 
evaluated  at  1 and  2 meters  for  each  wavelength  pair.  These  results 
are  shown  in  tables  5 and  6.  The  optimum  bands  for  this  case  appear  to 
be  .575-. 65  pm  or  .60-. 65  Pm.  These  are  very  close  to  the  actual  M-8 
bands  selected  for  processing  the  St.  Andrew  Bay  data  set  during  the 
1977  NCSC  experiment  [14],  The  performance  using  the  .60-. 65  pm  band 
pair  is  shown  in  Figure  7 for  sand  only  (upper  curve) ,and  averaged  over 
all  three  bottom  types  (lower  curve).  The  average  probability  of  correct 
classification  is  smaller  in  this  case  than  for  water  type  5 (Figure  6) 
because  of  the  similarity  of  the  reflectance  spectra  of  shoal  grass  and 
turtle  grass,  and  because  of  the  greater  attenuation  of  the  water.  When 
only  a discrimination  between  sand  and  vegetation  is  attempted,  the 
probability  of  correct  classification  is  much  higher,  as  shown  in  the 
upper  curve  of  Figure  7. 
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TABLE  3 

AVERAGE  CLASSIFICATION  ACCURACY  FOR  SAND,  MUD,  AND 
VEGETATION  IN  WATER  TYPE  5 AT  2 METERS  DEPTH 
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FIGURE  6.  AVERAGE  PROBABILITY  OF  CORRECT  CLASSIFICATION  OF  SAND.  MUD, 

AND  VEGETATION  IN  WATER  TYPES  3 AND  5,  FOR  BOTTOM  RECOGNITION 
ALGORITHM  TRAINED  IN  WATER  TYPE  5 AT  2 METERS  DEPTH. 
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AVERAGE  CLASSIFICATION  ACCURACY  FOR  SAND,  SHOAL  GRASS,  AND 
TURTLE  GRASS  IN  ST.  ANDREW  BAY  AT  1 METER  DEPTH 
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TABLE  6 

AVERAGE  CLASSIFICATION  ACCURACY  FOR  SAND,  SHOAL  GRASS, 
TURTLE  GRASS  IN  ST.  ANDREW  BAY  AT  2 METERS  DEPTH 
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FIGURE  7.  AVERAGE  PROBABILITY  OF  CORRECT  CLASSIFICATION  OF  SAND,  SHOAL 
GRASS,  AND  TURTLE  GRASS  (LOWER  CURVE),  AND  PROBABILITY  OF 
CORRECT  CLASSIFICATION  OF  SAND  ONLY  (UPPER  CURVE)  FOR  ST.  AND- 
REW BAY,  FLORIDA.  TRAINING  DEPTH  « 1 METER. 
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3.4  EMPIRICAL  EVALUATION 


The  two-channel  bottom  recognition  algorithm  has  been  applied  to 
two  actual  data  sets.  The  first  application  was  to  the  data  set  collected 
by  the  M-8  scanner  during  the  ERIM/NCSC  experiment  on  May  26,  1977  over 
St.  Andrew  Bay,  Florida  [14].*  Each  band  pair  in  this  data  set  was 
examined  for  maximum  separability  of  the  three  bottom  types  in  the  scene 
(sand,  shoal,  grass,  and  turtle  grass).  The  band  pairs  .545um  - .66um, 
,575um  - .66um,  and  ,61um  - ,66um  were  all  found  to  be  near-optimal  in 
terms  of  separation  distance.  Since  the  ,545pm  band  had  the  best  data 
quality  in  terms  of  noise  and  dynamic  range,  the  545  pm-. 66  pm  band  pair 
was  selected  for  processing.  The  transformed  signal  values  in  these 
bands  for  the  three  bottom  types  are  shown  in  Figure  8.  On  the  basis 
of  this  plot,  the  decision  boundary  separating  sand  and  shoal  grass  was 
taken  as  Y^  » 0.0,  and  the  decision  boundary  for  shoal  grass  and  turtle 
grass  was  taken  as  Y^  = -0.4. 

The  scene  shown  in  Figure  9 was  classified  using  this  set  of  para- 
meters. In  this  figure,  the  darkest  symbols  represent  pixels  classified 
as  turtle  grass,  the  intermediate  symbols  are  shoal  grass,  and  the  lightest 
symbols  are  sand.  The  blank  area  at  the  bottom  is  the  shoreline  near 
the  NCSC  Marina.  The  linear  feature  in  the  center  of  the  scene  is  a 
channel  having  a depth  of  4-5  meters.  The  area  surrounding  the  channel 
has  a depth  of  1-2  meters,  dropping  off  rapidly  to  5-6  meters  at  the 
top  of  the  scene.  The  dimensions  of  the  scene  are  approximately  320  meters 
in  each  direction  (there  is  some  distortion  of  scale  in  the  line  printer 
display) . An  aerial  photograph  of  the  scene  is  shown  in  Figure  10  for 
comparison. 

The  sandy  areas  near  the  shoreline,  in  scattered  patches,  and  on  the 
steep  slopes  of  the  channel  (center  of  scene)  and  at  the  edge  of  the  bank 
(top  of  scene)  appear  to  be  correctly  classified.  The  zones  of  shoal 
grass  in  the  center  and  of  turtle  grass  near  the  top  of  the  bank  appear 

* ' 

Data  collection  and  processing  for  this  experiment  were  carried  out  under 
contract  with  NCSC,  contract  no.  N61339-77-C-0059. 
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also  to  be  generally  correct.  The  classification  as  shoal  grass  in  the 
center  of  the  channel  and  in  deep  water  at  the  top  of  the  scene  is  more 
doubtful.  A complete  evaluation,  utilizing  subsurface  observations 
made  by  NCSC  personnel  at  the  time  of  the  overflight,  would  be  desirable 
(see  recommendations  for  further  study,  p.  6). 

The  second  application  of  the  bottom  recognition  algorithm  was  to 
Landsat  data  over  the  Bahamas.  A high  gain  Landsat  frame  (5249-14435) 
over  the  Great  Bahama  Bank  was  selected  for  analysis.  This  data  set 
had  been  previously  used  for  bathymetric  studies  in  a contract  with  the 
Defense  Mapping  Agency, *and  covers  an  area  which  is  currently  under 
intensive  study  as  a test  range  for  development  of  satellite  bathymetry 
techniques.  One  especially  interesting  part  of  the  scene  is  a long 
linear  feature  just  north  of  Mackie  shoal,  at  about  25'40'N  and  78°40'W. 
This  previously  uncharted  feature  was  visited  during  a field  trip  in 
October  1977,  and  was  found  to  consist  of  a sharp  depth  discontinuity 
of  about  2 meters,  with  a heavy  growth  of  turtle  grass  on  the  deeper 
(west)  side  and  sand  on  the  shallow  (east)  side.  A plot  of  Landsat  data 
values  in  MSS-4  (.50-.60ym)  and  MSS-5  (.60-. 70pm)  over  both  bottom  types  is 
shown  in  Figure  11.  Assuming  mean  oceanic  water  attenuation  coefficients, 
the  decision  boundary  between  sand  and  turtle  grass  is  about  Y^  * 3.0. 

A portion  of  this  linear  feature,  including  the  southern  terminus, 
is  shown  in  Figure  12.  The  dimensions  of  this  scene  are  about  6.9  km 
in  width  (left  to  right)  by  8.3  km  in  length  (top  to  bottom).  This  area 
was  processed  with  the  bottom  recognition  algorithm  using  mean  oceanic 
water  parameters.  Darker  symbols  represent  lower  values  of  and 
apparently  correspond  with  heavier  growths  of  turtle  grass.  Lighter 
symbols  represent  higher  values  of  Y^  and  apparently  correspond  with 
sandy  areas  in  the  scene.  Further  analysis  and  evaluation  of  this  scene 
is  recommended  to  determine  the  accuracy  and  depth  limits  of  the 
bottom  recognition  algorithm  as  applied  to  Landsat  data. 


*Contract  number  DMA-800-76-C-0057. 
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APPENDIX 

Passive  remote  sensing  techniques  for  mapping  water  depth 
and  bottom  features 


David  R.  Lyzertga 


Katin  processing  mrthiMls  are  reviewed,  anil  a new  method  is  proposed  for  ext  rai  l mg  water  depth  and  bot- 
tom  type  information  from  paaaive  nuiltispeetral  scanner  data  Limitations  of  each  technique  are  discussed, 
and  an  error  analysis  is  |>erformed  Using  an  analytical  model  for  the  radiance  over  shallow  water. 


Introduction 

Aerial  photography  of  shallow  water  areas  can  pro- 
vide useful  qualitative  information  on  bottom  compo- 
sition, the  distribution  of  henthir  algal  or  coral  com- 
munities, and  water  depth.  However,  the  interpreta- 
tion of  this  photography  is  impeded  by  the  fact  that 
water  depth  variations  are  not  easily  distinguished  from 
bottom  color  differences.  Surface  reflection  effects  add 
another  element  of  confusion  to  the  interpretation  of 
the  photography.  The  use  of  digitally  recorded  mul- 
tispectral  scanner  data  permits  corrections  to  be  made 
for  surface  reflection  effects  and  also  allows  the  possi- 
bility of  automatic  recognition  of  bottom  features  and 
water  depth  using  radiometric  techniques.  Past  re 
search  efforts  have  resulted  in  the  development  of 
specific  techniques  for  each  of  these  applications.1-2 
The  purpose  of  this  paper  is  to  discuss  the  limitations 
of  these  techniques  and  to  present  a more  general  set  of 
algorithms  for  both  applications. 

Ratio  Algorithms 

The  techniques  described  in  Refs.  1 and  2 were  de- 
veloped on  the  basis  of  a simple  water  reflectance  model 
which  accounts  for  the  major  part  of  the  signal  received 
by  a multispectral  scanner  over  clear  shallow  water,  but 
neglects  the  effects  due  to  scattering  in  the  water  and 
internal  reflection  at  the  water  surface.  According  to 
this  model,  the  radiance  in  a given  wavelength  band  (i) 
can  he  written  as 


L,  “ 4-  *,r„  exp(-»,/r),  to 

where  is  the  radiance  observed  over  deep  water  (due 
to  external  reflection  from  the  water  surface  and  scat- 
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tering  in  the  atmosphere);  k,  is  a constant  which  in- 
cludes the  solar  irradiance,  the  transmittance  of  the 
atmosphere  and  the  water  surface,  and  the  reduction 
of  the  radiance  due  to  refraction  at  the  water  surface; 
rp,  is  the  bottom  reflectance;  x,  is  the  effective  attenu- 
ation coefficient  of  the  water;  / is  a geometric  factor  to 
account  for  the  pathlength  through  the  water;  and  : is 
the  water  depth. 

The  simplest  method  of  extracting  water  depth  in- 
formation from  multispectral  scanner  data  is  to  invert 
Kq.  (1)  for  a single  wavelength  band.  An  extension  of 
this  method  would  lie  to  calculate  the  depth  from  two 
or  more  bands  and  average  the  results.  The  difficulty 
with  this  technique  is.  of  course,  that  changes  in  the 
bottom  reflectance  or  water  attenuation  cause  errors  in 
the  depth  calculation. 

The  water  depth  algorithm  developed  by  Polryn  ct 
al  1 relies  on  the  assumption  that  a pair  of  wavelength 
bands  can  be  found  such  that  the  ratio  of  the  bottom 
reflectances  in  these  two  bands  is  the  same  for  all  the 
hottom  types  within  a given  scene.  That  is.  for  bottom 
types  A.  R 


r*i  r, i3 


where  rA , is  the  reflectance  for  hottom  type  .4  in  hand 
1 etc.  The  water  depth  can  then  lie  calculated  from  the 
equation 


I 


to) 


where  R is  the  ratio  of  the  hot  tom -reflected  signals  in 
t he  two  hands: 


K ■ t/-i  — /.„»/(/.,  — l.,il  141 

If  the  assumption  expressed  by  Kq.  (2)  is  correct,  the 
depth  calculated  by  this  method  is  not  affected  by- 
changes  in  hottom  composition  in  the  scene.  The  depth 
is  also  insensitivp  to  changes  in  water  quality  if  the 
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difference  between  the  attenuation  coefficients  (»i  — *2) 
remains  constant.  In  many  cases  a pair  of  wavelengths 
can  be  found  for  which  Eq.  (2)  is  approximately  satis- 
fied, or  for  which  (n\  — *2)  remains  relatively  constant. 
However,  the  wavelengths  which  satisfy  one  criterion 
are  in  general  not  the  same  as  those  which  satisfy  the 
other,  and  if  changes  in  bottom  composition  or  water 
quality  are  too  large,  a pair  of  wavelengths  may  not  exist 
which  satisfies  either  criterion.  Nevertheless,  this 
method  has  been  used  with  some  success  for  extracting 
water  depths  from  both  satellite  and  aircraft  multis- 
pectral  scanner  data  over  relatively  clear  water  to  a 
depth  of  approximately  one  attenuation  length.1-4 

The  inverse  problem  is  to  extract  information  about 
the  bottom  reflectance  (or  bottom  composition)  from 
the  radiance  measured  by  the  multispectral  scanner. 
An  algorithm  was  developed  for  this  purpose'12  by  noting 
that  the  radiance  ratio  R should  be  independent  of  the 
water  depth  if  the  effective  water  attenuation  coeffi- 
cients are  the  same  in  both  bands.  From  the  simple 
water  reflectance  model  described  above,  this  ratio  then 
reduces  to 

ft  “ (*irsi)/(*jr«jt.  (5) 

which  may  be  used  as  an  index  of  the  bottom  type, 
provided  that  the  bottom  types  to  be  mapped  have 
different  reflectance  ratios  in  the  wavelength  bands 
selected. 

This  method  was  successfully  used  for  mapping  the 
distribution  of  Cladophora  (a  green  benthic  algae) 
under  a variable  depth  of  water  along  the  Lake  Ontario 
shoreline.2  The  primary  reason  for  the  success  of  this 
application  was  the  fact  that  the  vegetation  reflectance 
has  distinctive  features  (due  to  chlorophyll  absorption) 
in  the  blue-green  region  of  the  spectrum  where  the 
water  attenuation  is  at  a minimum.  Thus  it  was  pos- 
sible to  choose  two  bands  with  equal  water  attenuation 
coefficients  and  different  bottom  reflectance  ratios  for 
the  green  vegetation  and  the  sand  background. 

In  attempting  to  apply  this  method  to  a larger  num- 
ber of  water  types  and  bottom  materials,  however, 


X,  (0.47J  mkiI 

fit-  1-  Plot  of  ,Y  1 v»  .Y 2 for  water  type  .1  with  three  bottom  type*. 
Water  depth  rang**  from  0 m to  5 m on  each  curve. 


several  difficulties  arose.  First,  the  requirement  of 
equal  water  attenuation  coefficients  in  the  operating 
wavelength  bands  causes  operational  difficulties,  since 
the  band  positions  must  be  changed  when  operating  in 
different  water  types.  Second,  the  requirement  of 
different  bottom  reflectance  ratios  restricts  the  number 
of  bottom  materials  that  can  be  recognized,  since  ma- 
terials with  similarly  shaped  reflectance  spectra  (such 
as  sand  and  mud)  have  nearly  equal  reflectance  ratios. 
Finally,  both  ratio  methods  for  bottom  features  and 
water  depth  are  inherently  restricted  to  two  operating 
wavelength  bands.  Since  independent  information 
relating  to  both  water  depth  and  bottom  composition 
can  be  collected  simultaneously  in  several  wavelength 
bands,  methods  which  use  only  two  of  these  bands  do 
not  make  full  use  of  the  available  information. 


Mora  Ganaral  Algorithms 

In  order  to  develop  a more  general  set  of  algorithms 
for  water  depth  and  bottom  features,  the  simple  ra- 
diance model  described  above  was  modified  to  include 
the  effects  of  scattering  in  the  water  and  internal  re- 
flection at  the  water  surface  (see  Appendix  A).  Ex- 
amination of  this  model  shows  that  the  scattering  term 
has  the  same  depth  dependence  as  the  bottom-reflected 
radiance.  Except  for  the  effects  of  internal  reflect  ion, 
therefore,  the  total  radiance  for  the  case  of  direct  inci- 
dent radiation  may  be  written  in  the  same  form  as  Eq. 
(1),  with  the  actual  bottom  reflectance  replaced  by  an 
apparent  bottom  reflectance 


The  variables  on  the  right-hand  side  of  this  equation 
are  defined  in  Appendix  A.  This  exponential  depth 
dependence  suggests  the  use  of  the  transformation 

X,  ■ Intti  - L„).  (7) 

where  L„  is  the  deep-water  radiance  (including  scat- 
tering). The  purpose  of  this  transformation  is  to  li- 
nearize approximately  the  relationship  between  the 
transformed  radiances  and  the  water  depth.  The 
nonlinearities  which  remain  are  caused  by  internal  re- 
flection effects,  which  are  significant  only  for  very 
shallow  water  and  high  bottom  reflectances. 

A plot  of  the  transformed  radiances  over  three  bottom 
types  at  wavelengths  of  0.475  ym  and  0.525  *tm  is  shown 
in  Fig.  1.  The  water  parameters  used  in  this  example 
are  a composite  of  .lerlovV’  irradiance  attenuation 
coefficients  for  coastal  water  type  3 and  Petzold’s6 
scattering  parameters  for  Station  5 off  the  coast  of 
southern  California.  The  irradiance  attenuation 
coefficients  were  0.220  m-’  and  0.198  m_1,  and  the  total 
scattering  coefficient  was  0.275  m_l  at  the  wavelengths 
used.  In  the  absence  of  reliable  reflectance  measure- 
ments for  actual  bottom  materials,  the  reflectances  of 
beach  sand,  dark  soil  (representing  mud),  and  wheat 
leaves  (representing  aquatic  vegetation)  were  used  for 
this  calculation.  These  reflectances  are  plotted  in  Fig. 
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Fig.  2.  Spectral  reflectance*  of  aand,  mud.  and  green  vegetation. 


2.  For  these  water  and  bottom  parameters,  the  ra- 
diance due  to  scattering  is  equal  to  the  bottom-reflected 
radiance  at  a depth  of  about  6 m. 

Figure  1 represents  a simulated  data  set  for  an  ideal 
measurement  situation  with  no  noise  and  no  variation 
in  water  or  bottom  parameters.  For  an  actual  data  set, 
the  data  points  would  be  randomly  distributed  about 
the  lines  in  Fig.  1.  Except  for  the  slight  curvature 
caused  by  internal  reflection  effects,  the  transformed 
radiance  values  over  different  bottom  types  describe 
parallel  lines  in  X, -space  as  the  water  depth  varies 
continuously. 

For  an  N-band  system,  the  transformed  radiances  fall 
along  a set  of  parallel  lines  in  N-space.  Thus,  a second 
set  of  variables 

V,  - £ At]X,  (8) 

>•1 

can  be  obtained  by  rotating  the  coordinate  system  so 
that  the  Yjv  axis  is  parallel  to  this  direction  (see  Ap- 
pendix B).  If  the  linear  transformation  (8)  is  a pure 
rotation,  only  Ys  will  be  dependent  on  the  water  depth, 
while  all  the  other  variables  are  functions  only  of  the 
bottom  reflectance.  For  an  N-band  system,  this  results 
in  a set  of  N - 1 depth-invariant  signals  which  can  be 
used  as  inputs  to  a conventional  maximum  likelihood 
classification  algorithm.  The  remaining  variable  can 
be  written  as 

YN-Bm-  Ci.  (9) 

where  Bm  is  a function  of  the  bottom  composition  and 
C depends  only  on  the  water  attenuation  coefficients. 
If  the  bottom  material  can  be  recognized  by  the  proce- 
dure outlined  above,  and  the  value  of  Bm  determined 
for  each  bottom  type,  Eq.  (9)  can  then  be  used  to  cal- 
culate the  water  depth. 

Although  this  analysis  bears  a superficial  resemblance 
to  principal  component  analysis,7  the  only  actual  sim- 
ilarity is  the  use  of  a rotation  transformation.  In  the 
case  of  principal  component  analysis,  the  first  coordi- 
nate axis  is  aligned  in  the  direction  of  maximum  sample 
variance,  and  the  remaining  axes  are  aligned  in  the  or- 
thogonal directions  with  decreasing  sample  variance. 
The  purpose  of  the  principal  component  analysis  as 


applied  to  water  color*  is  to  reduce  the  number  of  vari- 
ables, since  only  a few  of  the  principal  components  are 
usually  needed  to  account  for  most  of  the  sample  vari- 
ance. In  the  analysis  presented  here,  the  directions  of 
the  transformed  coordinate  axes  are  not  necessarily 
related  to  the  sample  variance,  and  the  purpose  of  the 
analysis  is  not  to  reduce  the  number  of  variables  but  to 
remove  the  depth  dependence  from  all  but  one  vari- 
able. 

For  a t -band  system,  the  decision  rule  for  bottom 
classification  reduces  to  a simple  test  on  the  value  of  Y\. 
That  is,  the  bottom  is  classified  as  material  m if 

y».  < >'■  < yw2.  no) 

where  Ym|  and  Ym 2 are  the  lower  and  upper  decision 
boundaries,  respectively,  for  material  m.  If  the  radi- 
ances over  a given  bottom  type  m are  assumed  to  be 
normally  distributed  (due  to  system  noise)  with  mean 
value  Llm  and  standard  deviation  (NEAL),  in  band  i, 
the  probability  that  this  material  will  be  classified  as 
material  n is  approximately 

Pim.n)  « % erf  - 'h  erf  (^|)’  (ID 

where 

d,  - - >'„,)/»„.  (12) 

d*-  (?-  - y.j)/o„.  (13) 

- £ A 1,  (Nfc'Af.,)2/(£,  - L..)2  (141 
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The  total  probability  of  misclassification  for  material 
m is 

P(m)  » 1 - /’(m.m)  (15) 

The  total  probabilities  of  misclassification  for  the  three 
materials  in  Fig.  1 are  plotted  vs  depth  in  Fig.  3 for 
NEAL  = 0.05  mW  cm-2  sr-1  urn"1. 

It  can  be  shown  that  for  the  case  of  two  bands  with 
equal  attenuation  coefficients,  the  method  described 
above  is  equivalent  to  the  ratio  method.  However,  in 
this  case  the  performance  of  both  methods  is  worse  than 
that  shown  in  Fig.  3 because  of  the  similar  reflectance 


WATER  DEPTH  Imwer*) 

Fig.  3.  Total  probability  of  misclassification  for  three  bottom  type* 
in  water  type  3,  using  the  proposed  method. 
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ratios  of  sand  and  mud.  Thus,  the  performance  of  the 
ratio  method  can  always  lie  equaled  and  can  usually  be 
ex<-cedcd  by  the  use  of  the  proposed  method  with  a more 
optimum  pair  of  wavelengths.  Further  improvement 
in  bottom  classification  would  be  expected  if  more  than 
two  (>ands  were  used. 

After  the  bottom  classification  has  been  performed, 
the  water  depth  can  lie  calculated  from  Kq.  (9),  using 
the  appropriate  value  of  H„,  for  material  m.  Depth 
errors  arise  from  the  variance  of  Y\  due  to  noise,  which 
is  approximately  given  by 

U >'*>*•  1 )'•/{[.,  -/.„)*  (16) 

and  from  errors  in  bottom  classification.  Thus  the  total 
rms  depth  error  for  Isittom  type  m is 


UV.vl2  + /’Im.n  ){Hm  - 8„  1- 


The  average  depth  error  for  the  three  bottom  types 
in  water  type  3 is  plotted  vs  depth  in  Fig.  4.  For  com- 
parison. the  depth  error  due  to  noise  only  using  the  ratio 
method  is  also  plotted  in  Fig.  4.  The  rms  depth  error 
due  to  noise  for  the  ratio  method  is  given  by 

a/,  • — - — I — Y ♦ ( y |‘  m 

The  depth  error  is  in  general  minimized  by  choosing 
wavelength  hands  with  the  smallest  attenuation.  In  thi 
ratio  method,  however,  the  sensitivity  to  noise  increases 
rapidly  as  the  difference  between  the  attenuation 
coefficients  in  the  two  bands  decreases.  Since  this 
difference  is  only  0.025  m 1 in  the  above  example,  these 
are  not  the  optimum  wavelengths  for  the  ratio  method. 
For  any  choice  of  wavelengths,  however,  the  depth  error 
due  to  noise  is  larger  for  the  ratio  method  than  for  the 
method  described  above.  The  ratio  method  is  also 
subject  to  errors  due  to  changes  in  bottom  composition: 
although  these  can  lie  reducer!  by  an  appropriate  choice 
of  wavelengths,  they  can  never  in  practice  be  completely 
eliminated. 


Conclusions 

The  ratio  algorithms  for  water  depth  and  bottom 
features  mapping  are  relatively  simple  techniques  which 
give  acceptable  results  in  many  situations.  However, 
these  algorithms  are  limited  by  operational  restrictions 
which  reduce  their  applicability  and  utility.  A more 
general  algorithm  has  been  defined,  and  a preliminary 
evaluation  of  this  method  has  been  performed  using  a 
simulation  model  for  the  water  radiance. 

The  advantages  of  this  method  for  mapping  bottom 
features  are  (1)  increased  operational  flexibility,  in  that 
the  wavelength  bands  are  not  limited  to  those  with  equal 
water  attenuation  coefficients,  (2)  better  discrimination 
of  bottom  materials  with  similarly  shaped  reflectance 
spectra,  and  (3)  improved  performance  through  the  use 
of  more  than  two  wavelengths  bands. 

The  advantages  of  this  method  for  mapping  water 
depth  are  ( 1 ) increased  operational  flexibility,  since  the 
wavelength  bands  are  not  restricted  by  the  requirement 
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KiC  1 Total  depth  error  using  the  proponed  method  and  depth  error 
due  to  ihhm*  using  the  ratio  method  for  water  and  Imttom  types  shown 
in  Fig  I 


of  equal  bottom  reflectance  ratios  for  all  bottom  ty|>e.s, 
(2)  lower  sensitivity  to  noise,  and  (31  improved  perfor- 
mance through  the  use  ot  more  than  two  wavelength 
bands. 

The  disadvantage  of  the  algorithm  defined  here  is 
that  it  is  more  complex  and  therefore  somewhat  more 
difficult  to  implement  than  the  ratio  me  hods.  The 
subtraction  and  division  operations  required  for  the 
ratio  methods  can  be  implemented  by  either  analog  or 
digital  processors,  whereas  the  method  defined  here 
requires  digital  computation.  The  input  parameters 
for  this  algorithm  are  also  somewhat  more  difficult  to 
determine  than  those  for  the  ratio  methods. 

Numerical  results  illustrating  the  performance  of  this 
method  have  tieen  presented  for  one  example  situation. 
^ hese  results  should  not  be  considered  as  definitive  of 
the  lwst  possible  performance,  since  the  wavelengths 
considered  are  not  necessarily  the  optimum  ones  nor  are 
the  reflectances  used  necessarily  representative  of  ac- 
tual bottom  types.  In  addition,  lie! ter  performance  can 
lie  obtained  by  reducing  system  noise  and  by  using  more 
than  two  wavelength  bands.  The  value  of  .Vf.’A/.  used 
in  these  calculations  was  0.05  mW  cm'2  sr'1  pin'1. 
This  is  a typical  value  for  an  aircraft  multispectral 
scanner  with  an  angular  field  of  view  of  a few  millirad- 
ians,  a spectral  resolution  of  about  0.025  pm,  and  a t ol- 
lector  area  of  about  80  cm2.  This  noise  equivalent  ra 
diance  can  readily  lie  reduced  by  a factor  of  2 or  3,  with 
a corresponding  reduction  in  the  depth  error,  by  in- 
creasing the  collector  area  or  reducing  the  spatial  or 
spectral  resolution  of  the  system. 

This  work  was  supported  by  the  Office  of  Naval  Re- 
search, contract  N000- 14-74-0-0273. 

Appendix  A:  Shallow-Water  Radiance  Model 

The  radiances  shown  in  Fig.  I were  calculated  using 
a combined  water  atmosphere  radiance  model  which 
includes  the  effects  of  scattering  in  the  atmosphere  and 
reflection  at  the  water  surface,  as  well  as  the  compo- 
nents originating  in  the  water  itself.  The  atmospheric 
effects  are  calculated  from  the  double-delta  approxi- 
mation,9 using  the  atmospheric  parameters  tabulated 
by  Klterman.10  A solar  zenith  angle  of  20°  and  a plat 
form  altitude  of  1 km  were  used  for  the  calculations  in 
this  report. 
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The  water  radiance  is  calculated  using  the  quasi- 
single scattering  approximation,11  modified  to  include 
the  effects  of  reflection  from  the  bottom12  and  internal 
reflect  ion  at  the  water  surface. 1:1  The  effects  of  internal 
reflection  are  included  to  all  orders  by  assuming  that  the 
upwelling  radiance  distribution  just  lieneath  the  water 
surface  has  a uniform  angular  distribution.  Thus,  the 
upwelling  underwater  radiance  can  be  written  as 

Eo'Siu'.uo.t)  + f 
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Jo  Jo 
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where  Eo  is  the  direct  solar  irradiance  penetrating  the 
water  surface,  //(/i',0")  is  the  transmitted  skv  radiance 
below  the  surface,  no’  is  the  cosine  of  the  apparent  solar 
zenith  angle  below  the  surface,  R(n")  is  the  Fresnel  re- 
flectance of  the  water  surface,  and 

Nt».v'.iM  - j;  U + s')-1  1 1 - exp  | Kz  j | d(x.) 

+7exp["(;+^)#f*]-  ,A2‘ 

where  K is  the  irradiance  attenuation  coefficient,14  r 
is  the  water  depth,  dip., ) is  the  volume  scattering  func- 
tion, rn  is  the  Ixittom  reflectance,  and 

(i,  = -pp'  + |(l  - p2)(l  - m'2)!1'2  cos «.  (A3) 

'I'he  results  of  this  model  agree  with  exact  calculations 
using  the  Monte  Carlo  method1''  to  within  about  10% 
tor  values  of  the  single-scattering  albedo  less  than  0.8 
and  bottom  reflectances  less  than  50%. 


Appendix  B:  Coordinate  System  Rotation  Parameters 

The  assumption  involved  in  making  the  coordinate 
transformation  (8)  is  that  the  variables  A’,  are  linearly 
correlated  with  the  water  depth  z,  that  is, 

X,  = a,  - b ,z  (Bl) 

with  possibly  different  values  of  a,  for  each  Ixittom  type. 
It  is  desired  to  transform  this  set  of  N depth-dependent 
variables  into  a new  set  of  N - 1 deptb-independent 
variables  (V|  ..  V'/v- 1 ) and  one  depth-dependent 

variable  (Vjv).  The  required  transformation  may  be 
visualized  as  a pure  rotation  in  N -space  such  that  the 
Y'v  axis  is  oriented  in  the  direction  specified  by  the 
parametric  Eq.  (Bl).  Thus,  the  transformation  pa- 
r imeters  for  V.\  are 


The  remaining  parameters  can  be  calculated  using 
the  conditions  for  the  orthonormality  of  the  new  coor- 
dinate axes 


* 10,  i * J 

1 AtkA,k  * . (B 3) 

km  1 l 1,  I * J 

and  the  condition  for  a proper  rotation 

(B4I 

These  conditions  uniquely  define  the  transformation 
matrix  only  for  the  case  N = 2.  For  three  or  more  di- 
mensions, however,  a unique  solution  can  be  obtained 
by  requiring  that  the  coefficients  for  Y\. . .Yn-  i be  the 
same  for  the  (N  + 1 ^dimensional  case  as  for  the  N- 
dimensiona!  case.  The  coefficients  for  the  (N  — 1) 
depth-independent  variables  are  then 

A0  * 6,*i6,  ^ £ bk2j  bk  -j  for  ;*  1 . . . (B5) 

At)  * - ^ V ^ 6a*  j for  ;■  i + 1,  (B6) 

Al}  * 0 for;  *i  + 2 N.  (B7) 

These  equations,  along  with  Eq.  (B2),  completely 
define  the  required  coordinate  transformation.  The 
values  of  h,  can  be  calculated  if  the  water  attenuation 
coefficients  are  known  or  can  be  empirically  obtained 
from  a regression  analysis  of  the  measured  radiance 
values  over  a uniform  bottom. 
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